arXiv: 1508.00101 vl [gr-qc] 1 Aug 2015 

O.M n CJ Ci-K-(fa M-srl-r-r'n. “* ^ 


Application of Sturm’s theorem to marginal stable circular orbits 
of a test body in spherically symmetric and static spacetimes 


Toshiaki Ono^, Tomohito Suzuki^, Naomasa Fushimi^, Kei Yamada^ and Hideki Asada^ 

^ Faculty of Science and Technology, Hirosaki University, Hirosaki 036-8561, Japan 
^ Department of Physics, Kyoto University, Kyoto 606-8502, Japan 


PACS 04.25.Nx - Post-Newtonian approximation; perturbation theory; related approximations 
PACS 95.30. Sf - Relativity and gravitation 
PACS 04.70. -s - Physics of black holes 

Abstract - In terms of Sturm’s theorem, we reexamine a marginal stable circular orbit (MSCO) 
such as the innermost stable circular orbit (ISCO) of a timelike geodesic in any spherically sym¬ 
metric and static spacetime. MSCOs for some of exact solutions to the Einstein’s equation are 
discussed. Strum’s theorem is explicitly applied to the Kottler (often called Schwarzschild-de 
Sitter) spacetime. Moreover, we analyze MSCOs for a spherically symmetric, static and vacuum 
solution in Weyl conformal gravity. 


Introduction. — In general relativity, the orbital ra- 
lius has a lower bound that is called the innermost sta- 
)le circular orbit (ISCO). The existence of the ISCO in 
he Schwarzschild solution for the Einstein’s equation is 
ascinating [I]. Moreover, ISCOs may play key roles in 
Lstrophysics as well as in gravity theory. For instance, 
SCOs are of great importance in gravitational waves as- 
ronomy [5], because ISCOs are thought to be the location 
d the transition from the inspiralling phase to the merging 
)ne, especially when a compact object is orbiting around 
I massive black hole probably located at a galactic cen- 
er. Furthermore, in high energy astrophysics, ISCOs are 
elated to the existence for the inner edge of an accretion 
lisk around a black hole [3]. It is thus expected that a 
measurement of the ISCO radius will bring us important 
information on the strong gravity, especially the nonlinear 
spacetime geometry that is beyond the solar-system tests. 
For near-future astrophysical tests of the no-hair theorem 
for black holes, therefore, it is intriguing to study the na¬ 
ture of the ISCOs in non-Schwarzschild spacetimes such 
as black holes with electric charges and/or scalar fields [1] 
and black holes in modified gravity theories [5] . 

It has been thought that unstable circular orbits would 
be generally produced by strong gravity in general rela¬ 
tivity. This is the case for the Schwarzschild spacetime. 
If the cosmological constant for instance is added into the 
Einstein’s equation, however, it is not always the case. 
Stuchlik and Hledik [5] have pointed out that the outer¬ 
most stable circular orbit (OSCO) of a test body is possi¬ 


ble in the Kottler (often called the Schwarzschild-de Sitter) 
spacetime [?]■ The ISCO and OSCO are a boundary be¬ 
tween a stable region and an unstable one. Hereafter, we 
call it a marginal stable circular orbit (MSCO). Note that 
the wording “marginally stable” has a different meaning 
in the theory of dynamical systems. Hence, this paper 
prefers to use the wording “marginal stable circular or¬ 
bits” . If MSCOs are two and only two in a spacetime, one 
MSCO may correspond to the ISCO and the other may 
correspond to the OSCO. Note that the number of MSCOs 
might be larger than two, if such a spacetime geometry is 
described by a very complicated form of the metric. For 
this case, the smallest MSCO may be called the ISCO and 
the largest one may be the OSCO. 

Even counting the number of MSCOs for a non- 
Schwarzschild metric, though it is apparently elementary, 
is not a straightforward task. The purpose of this brief 
paper is to reexamine, in terms of Strum’s theorem [5], 
MSCOs of a timelike geodesic in any spherically sym¬ 
metric and static spacetime that may have a deficit an¬ 
gle. MSCOs for some of exact solutions to the Einstein’s 
equation are mentioned: Schwarzschild, Kottler (often 
called Schwarzschild-de Sitter) [7], Reissner-Nordstrom 
(RN) i, Janis-Newman-Winicour (JNW) spacetimes [TU], 
and Barriola-Vilenkin monopole with mass that is closely 
related to a static string m We illustrate how Strum’s 
theorem is explicitly applied to the Kottler spacetime. 

Moreover, we analyze a spherically symmetric, static 
and vacuum solution in Weyl conformal gravity, which is 
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an alternative gravity model. The exterior metric to a 
static spherically symmetric distribution in Weyl confor¬ 
mal gravity was found by Mannheim and Kazanas m- 
This solution expresses a black hole for some parameter 
range and it is a naked singularity for another parame¬ 
ter region. Recently, the Weyl-gravity solution has been 
studied to explain rotation curves of galaxy samples with¬ 
out introducing a dark matter component m and the 
model may provide also a possible size of individual galax¬ 
ies [14]. Moreover, the gravitational lensing by this black 
hole model has been discussed by several groups [TSlfTB] . 

Throughout this paper, we use the unit of G = c = 1. 

Application of Strum’s theorem to exact solu¬ 
tions for the Einstein’s equation. — 

Equation for a location of a MSCO. We consider 
spherically symmetric and static spacetimes. A general 
form of the line element for these spacetimes is 

ds^ = -A{r)dt^ + B{r)dr^ + C{r){d9^ + sin^ 0#^), (1) 

where we assume gtt = —A(r) < 0, grr = B{r) > 0, 
ggg = C(r) > 0, namely the spacetime signature as 
—, -b, -b, -b. One may often choose C{r) = r^. However, it 
would be more convenient to keep G(r) for some cases in 
the literatures where C (r) [TTUfTT] , because we do not 
need to add calculations of the coordinate transformation 
for reaching G(r) = for such cases. 

A necessary condition for the existence of a MSCO is 
known as [ISlfTO] 

dr \A(r )) dr^ \G(r) ) dr \G(r) ) dr'^ \^(''’) / 

( 2 ) 

This equation determines the radius of the MSCO, if the 
MSCO exists. Hereafter, we call Eq. (21) MSCO equation. 
Note that the MSCO equation does not contain a metric 
component B{r) nor the constants of motion as E and L. 

Given a root for Eq. E^ and are expressed as 


1 d 1 

A dr 


(3) 

1 d , 

( 1 A 

(4) 


VA(r)j’ 


where we define a determinant as 


1 1 



Any root r for Eq. (2|) can be substituted into Eqs. (O 
and in order to see whether it satisfies the sufficient 
condition as 0 < < 00 and 0 < < 00 . If the 

sufficient condition is satisfied, this r is a MSCO radius, 
denoted as cm SCO- If not, the root is unphysical and it 
must be discarded. Note that the expressions of E"^ and 
in Eqs. m and are applicable not only to MSCOs 
but also to any circular orbit. 


MSCOs in terms of Strum’s theorem . In this subsec¬ 
tion, we apply Strum’s theorem to some of exact solutions 
of the Einstein’s equation in order to classify the number 
of MSCOs. 

Schwarzschild spacetime: 

Let us begin with the Schwarzschild spacetime as 

ds^ = — fl — —) dt^-\ — „ -br^(d0^-bsin^ 9d(j)^), (6) 
V r J 1 - -A 
r 

where the Schwarzschild radius is defined as rg = 2M for 
the ADM mass M. Eor this spacetime metric, Eq. @ 
becomes linear as 

r-Mg = 0, (7) 

where we assume r 0 because r = 0 is the spacetime 
singularity. Hence, we obtain rMSCO = 3rg. 

Kottler (Schwarzschild-de Sitter) spacetime: 

Next, we consider the Kottler spacetime [Tj. The line ele¬ 
ment is 



-br^ {d9'^ -b sin^ 9d(lP ), (8) 

where A denotes the cosmological constant. Note that 
this spacetime is not asymptotically flat. The timelike 
geodesics in the Kottler spacetime have been often studied 
by several authors 20], especially on the periastron shift, 
though none of them has examined the allowed region of 
the orbits. 

Eor this spacetime, the MSCO equation becomes 

^Ar^ — 5rgAr^ — rgr -b = 0, (9) 

where we assume r 7^ 0 and r ^ ex. It is convenient to 
use the normalized variables in terms of the Schwarzschild 
radius rg. We define x = rjrg and A = Ar^/S. The above 
quartic equation becomes 

8Ax^ - I5Aa;3 - a;-b 3 = 0. (10) 

The remaining task is to investigate positive zeros for 
this quartic equation as a necessary condition for MSCOs. 
A general formula for solving a quartic equation such 
as Cardano’s one gives a complicated form of the root. 
Hence, it would be preferred to use an alternative method. 

Strum’s theorem [5] expresses the number of distinct real 
roots of a polynomial p located in an interval in terms of 
the number of changes of signs of the values of the Strum’s 
sequence (called also the Sturm’s chain) at the bounds of 
the interval. 

Sturm’s theorem: Let p{r) denote a polynomial. Apply¬ 
ing Euclid’s algorithm to p{r) and its derivative, Sturm’s 
sequence is constructed as 

Poic) = p{r), 
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PiW = P'o{r), 

P2{r) = Pi{,r)qQ{r) - po{r), 

Psir) = P2{r)qi{r) - pi{r), 

0 = Pn{r)q„-iir) - pn-i{r), (11) 

where qi{r) is the quotient of Pi{r) by pi+i{r). Let V(a) 
denote the number of the sign changing (ignoring zeros) in and 
the Sturm’s sequence at r = a. Then, V{a) — V{b) gives the 
number of distinct roots of p{r) between a and b, where 
a < b. 

Strum’s sequence for Eq. (na in the Kottler spacetime 


is obtained 

as 


Mr) = 

SAr"^ — 15rgAr^ — r^r + 3rg, 

( 12 ) 

Mr) = 

32Ar^ — 45rgAr^ — r^. 

(13) 

Mr) = 

^rl{225XrM32rgr-l23rl), 

128rl 

(14) 

Mr) = 

9 

50625A 



X [-2(128 + 4725A)r- + 3rg(328 + 3375A)], 

(15) 

1 t;i 

= - 512(128 

(16) 

where / 4 (r) is a constant for the present case. Henceforth, 
we focus on non-Schwarzschild case (A 0). 

First, we consider the case of A > 0, for which we have 

sgn[/o( 0 )] > 0 , 
sgn[/i( 0 )] < 0 , 
sgn[/ 2 ( 0 )] < 0 , 
sgn[/ 3 ( 0 )] > 0 , 

sgn[/ 4 ( 0 )] = sgn[-A(4 - 27A)(16 - 16875A)],(17) 

and 

sgn[/o(oo)] > 0 , 
sgn[/i(oo)] > 0 , 
sgn[/ 2 (oo)] > 0 , 
sgn[/ 3 (oo)] < 0 , 

sgn[/ 4 (oo)] = sgn[—A(4 — 27A)(16 — 16875A)],(18) 

where sgn denotes the sign function. 

(1) If 16/16875 < A < 4/27, the number of the sign 
changing is V(0) = 2 and V(oo) = 2, so that there cannot 
exist any positive root. (2) If 0 < A < 16/16875 or A > 
4/27, the sign changing number is V(0) = 3 and V(oo) = 
1 , which lead to two positive roots. (2a) If A > 4/27, 
however, < 0, which means no MSCO. (2b) On the 
other hand, both > 0 and > 0 are satisfied by 
0 < A < 16/16875, which means two MSCOs. 


Next, we consider A < 0. Then, we have 

sgn[/o(0)] > 0, 
sgn[/i(0)] < 0, 
sgn[/ 2 ( 0 )] < 0, 

sgn[/ 3 ( 0 )] =sgn[-(328 + 3375A)], 
sgn[/ 4 ( 0 )] > 0, (19) 

sgn[/o(oo)] < 0 , 
sgn[/i(oo)] < 0 , 
sgn[/ 2 (oo)] < 0 , 

sgn[/ 3 (oo)] = sgn[128 + 4725A], 
sgn[/ 4 (oo)] > 0 . ( 20 ) 

Regardless of sgn[/ 3 ( 0 )] and sgn[/ 3 (oo)], we find V(0) = 
2 and V(oo) = 1, which lead to a single positive root. One 
can also see that it satisfies > 0 and > 0. If A < 0, 
then, there always exists a single MSCO corresponding to 
the ISCO. 

Strum’s theorem (as the necessary condition) and the 
positive and (as the sufficient condition) tell that 
there are four cases: 

Case 1: A = 0. There is the single positive zero cor¬ 
responding to the ISCO. This is in agreement with the 
Schwarzschild case. 

Case 2: 0 < A < 16/16875. We have two positive zeros, 
namely two MSCOs, where one is corresponding to the 
ISCO and the other is the OSCO. 

Case 3: 16/16875 < A. There is no MSCO (after the ISCO 
and the OSCO merge at A = 16/16875). This implies that 
every circular orbit becomes unstable for this case. This 
number 16/16875 was found by Stuchlik and Hledik [5], 
where this is expressed as y = 12/15^ in Eq. (24) of their 
paper. 

Case 4: A < 0 (anti-de Sitter case). There is only the 
single MSCO. 

Other examples in general relativity: 

Let us briefly mention MSCOs in Reissner-Nordstrom 
(RN) [5], Janis-Newman-Winicour (JNW) spacetimes [TU] . 
and Barriola-Vilenkin monopole with mass m- MSCO 
equations are cubic, quadratic and linear, respectively, 
though at first glance the JNW case appears to lead to 
a higher order polynomial m- Therefore, Strum’s theo¬ 
rem is not needed for these cases. 

Spherically symmetric, static and vacuum solution 
in Weyl conformal gravity: 

Finally, we study a spherically symmetric, static and vac¬ 
uum solution in Weyl conformal gravity m as an alter¬ 
native gravity model. The exterior spacetime metric in 
vacuum regions was found by Mannheim and Kazanas m 
as 

= -A{r)dt^ + + sin^ 0d<f), 

A{r) = 1 — 3/3^ — — - — (21) 
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where /3, 7 and k are the integration constants to the vac¬ 
uum equation in Weyl conformal gravity m- The k pa¬ 
rameter may play the same role as the cosmological con¬ 
stant in the Kottler spacetime. 

It is more convenient to rearrange this original form of 
the metric into another one m as 

= -A(r)dt^ + -^dr^ + r^ldO^ + sin^ 0d(j)'^), 
Air) 

^TYl 

A{r) = a -hyr —fcr^, (22) 

r 

where we defined a = ^/1 — Gmy, and m = /3(2 — 3 / 37 )/ 2 . 
The negative mass case prohibits a circular orbit because 
of its repulsive force. This paper thus focuses on m > 0 
to study MSCOs. 

In this paper, both cases of 7 > 0 and 7 < 0 are con¬ 
sidered. When we wish to discuss a possible explanation 
of galactic rotation curves in Weyl gravity, 7 > 0 is pre¬ 
ferred [HITS]. On the other hand, 7 parameter cannot be 
constrained by galactic rotation curves, if one takes the 
standard cosmological point of view that galactic rotation 
curves can be explained by introducing a dark matter com¬ 
ponent. In this stance, therefore, the sign of 7 parameter 
can be uncertain at present. 

Let us begin with the most general case (fc 0 and 
7 ^ 0 ). Straightforward calculations show that MSCO 
equation becomes quintic in r, where the fifth-order term 
is due to a coupling between k and 7 parameters. The 
quintic equation is not algebraically tractable. Moreover, 
roughly speaking, terms that may play a role at cos¬ 
mological distances are negligible at smaller scales such as 
galactic and stellar sizes compared with the linear terms. 
In this paper, therefore, we focus on /c = 0 and 7 0 in 

order to neglect the terms in the metric. This case in 
Weyl conformal gravity may be used to discuss a modifi¬ 
cation to a stellar-mass black hole model, because k can 
be effective only at cosmological distances as mentioned 
above. See also the metric for Kottler solution. 

Interestingly, Eq. (EH) with k = 0 and 7 7 ^ 0 ac¬ 
cidentally coincides with an exact solution of Einstein 
equation with the quintessential matter surrounding a 
black hole (often called Kiselev black hole [H]), if the 
parameter in Kiselev black hole can be adjusted suitably 
(3w„-|-l = — 1). Note that the parameters of Weyl-gravity 
black holes and Kiselev ones have very different physical 
origins: The parameters in Weyl-gravity black holes are 
the integration constants coming from a higher derivative 
theory, while those in Kiselev black holes are due to the 
equation of state for a dark energy component, more pre¬ 
cisely the potential form that the quintessence field obeys. 

According to Eqs. dH and (0]), > 0 and > 0 for 

this spacetime metric are 

7 r^ -I- 2-\/l — 6mr — 6m > 0, (23) 


and 


-I- 2 m 

+ 2 1/1 — 6mr — 6m 


> 0 , 


(24) 


respectively. 

If 7 > 0, A{r) = 0 gives the black horizon radius as 


rn 


VT+'Smy — ^/T^^6mrf 

2 ^ 


If 7 < 0, A{r) = 0 gives the black horizon radius as 


(25) 


^ Hin 


^ Hout 


•v/W-'Gmy -|- y/1 + 2m'y 

2 ^ ’ 

y/1 — Gmy — + 2m"f 

2 ^ ’ 


(26) 

(27) 


where the former and the latter correspond to the inner 
horizon and the outer one, respectively. Therefore, the 
present metric is a black hole solution if and only if 



1 1 

— -— < 7 < -—. 

2 m 6 m 

(28) 

The MSCO 

equation for k = 0 and 7 7 ^^ 0 becomes 



7 ^r"^ -I- 37 -^/l — bmyr^ — 12m7r^ 



+2my/l — 6mjr — 12 m^ = 0 . 

(29) 

Let us define 




7 = 2 m 7 , 

(30) 



(31) 

so that the Sturm’s sequences can be in a simpler form as 

fo{r) 

= — 7 ^r"’’ — 37\/1 — 37r^ -|- 677 -^ 



-•\/l - 377 - 1 - 3 , 

(32) 

h(.r) 

= —47^7^ — 9-\/l — 3777 ^ + 1277 



-y'l - 37 , 

(33) 

f2{r) 

= --^(9-ll7)^^+3v^l-377 

16 



-t|^(1 + 137), 

(34) 

fsir) 

- (9-113)79 787 + 257^)r 




(35) 

U{r) 

3 ( 9 - 117)2 

167 ( 3 - 257)2(-3-f 7)2 



x(l-f7)(l-7 907 - 2372 ), 

(36) 


where / 4 (r) is a constant. Eq. (1^51) is rewritten as 


-1<7<3. 


(37) 


Case 1: 7 > 0. The sign of the Sturm’s sequence for this 
case is 


sgn[/o( 0 )] 

> 

0 , 

(38) 

sgn[/i( 0 )] 

< 

0 , 

(39) 

sgn[/ 2 ( 0 )] 

< 

0 , 

(40) 

sgn[/ 3 ( 0 )] 

< 

0 , 

(41) 

sgn[/4(0)] 

= 

sgn[l-7 907- 2372]. 

(42) 
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and 


sgn[/o( 

00 )] 

< 0, 

(43) 

sgn[/i( 

00 )] 

< 0, 

(44) 

sgn[/2( 

00 )] 

= sgn[-(9-ll7)]. 

(45) 

sgn[/3( 

00 )] 

= sgn[9-787-b257^], 

(46) 

sgn[/4( 

00 )] 

= sgn[l-b 907 - 237 ^]. 

(47) 


There is a constraint as 7 < 1/3, because there exists 
the — 37 term in the metric. Hence, we can show 

sgn[/ 2 (oo)] < 0, (48) 

sgn[/ 4 ( 0 )] = sgn[/ 4 (oo)] 

> 0. (49) 

Regardless of sgn[/ 3 (oo)], hence, V(0) = 2 and V(oo) = 
1. Therefore, there exists only the single positive root. 
Moreover, it satisfies > 0 and > 0. Hence, there is 
the only one MSCO, namely ISCO. 

Case 2: 7 < 0. This case is rather tricky. Non-trivial 
calculations are needed as shown below. The sign of the 
Sturm’s sequence for this case is 


sgn[/o(0)] 

> 

0, 

(50) 

sgn[/i(0)] 

< 

0, 

(51) 

sgn[/2(0)] 

= 

sgn[l -b 137], 

(52) 

sgn[/3(0)] 

> 

0, 

(53) 

sgn[/4(0)] 

= 

sgn[-(l -b 7)(1 -b 9O7 ' 

- 237 ')], (54) 


and 


sgn[/o(oo)] < 

0, 


(55) 

sgn[/i(oo)] < 

0, 


(56) 

sgn[/ 2 (oo)] < 

0, 


(57) 

sgn[/ 3 (oo)] = 

sgn[- 

-(9-787+ 257 ')], 

(58) 

sgn[/ 4 (oo)] = 

sgn[- 

-(l-b7)(l-b 907 - 2372 )]. 

(59) 


For — 1 < 7 < 0, one can show that 

sgn[-(l + 7 )] < 0 , (60) 

sgn[/ 3 (oo)] < 0 . (61) 

Regardless of the sign of / 2 ( 0 ), therefore, the number of 
the positive roots changes at 


As a result, there are three types for — 1 < 7 < 0: 

(2-1) If (45 — 32-\/2)/23 < 7 < 0, there are three positive 
zeros. However, one of them is located outside the outer 
horizon, while the others are placed between the two hori¬ 
zons and they satisfy both > 0 and > 0. Namely, 
there are two MSCOs. 

(2-2) If 7 = (45 - 32^2)723, the ISCO and the OSCO 
merge. 




r 


Fig. 1: 7 -dependence of the MSCO radius. The horizontal 
axis denotes the MSCO radius tmsco as a root for the MSCO 
equation. The vertical axis denotes 7 parameter. The shaded 
parts denote prohibited regions where A'^< 0 , I/^< 0 , 7<—1 
or 7 > 1/3. The dashed curve in the shaded regions denotes 
the horizon location for a given 7 value. Top: 7 € [—1.2, 0.5]. 
Bottom: 7 € [—0.1, 0] for a zoom-up figure of the Case 2. 


(2-3) If -1 < 7 < (45 - 32^2)723, the MSCO equation 
has one positive root. However, it does not satisfy E'^ > 0. 

Bringing the Case 1 (corresponding to 7 > 0) and the 
Case 2 (corresponding to 7 < 0) together. Figure [T] shows 
the dependence of the MSCO radius on 7 . 

Conclusion. — We reexamined, in terms of Strum’s 
theorem, MSCOs of a timelike geodesic in any spherically 
symmetric and static spacetime that may have a deficit an¬ 
gle. MSCOs for some of exact solutions to the Einstein’s 
equation were mentioned: Schwarzschild, Kottler (often 
called Schwarzschild-de Sitter) [ 7 ], Reissner-Nordstrom 
(RN) i, Janis-Newman-Winicour (JNW) spacetimes [TH] . 
and Barriola-Vilenkin monopole with mass We il¬ 

lustrated how Strum’s theorem is explicitly applied to 
the Kottler case. This suggests that Strum’s theorem is 
widely applicable for classifying MSCOs for some space- 
time, when the MSCO equation under study is equivalent 
to a polynomial. 
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Moreover, we analyzed MSCOs for a spherically sym¬ 
metric, static and vacuum black hole solution in Weyl 
conformal gravity, where we neglected the square terms by 
assuming fe = 0. MSCO equation becomes quintic for a 
more general black hole solution in Weyl gravity with non¬ 
zero k parameter. Numerical investigations are needed to 
classify MSCOs in this quintic case. It is left as a future 
work. 

* * * 
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